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ABSTRACT 

We have serendipitously detected a strong O Vl-bearing Lyman limit system at Zats = 0.3558 toward 
the QSO J1009+0713 (zem = 0.456) in our survey of low-redshift galaxy hales with the Hubble Space 
Telescope's Cosmic Origins Spectrograph. Its total rest-frame equivalent width of Wr = 835 ± 49 mA 
and column density of log A^(0 VI) = 15.0 are the highest for an intervening absorber yet detected in 
any low-redshift QSO sightline, with absorption spanning at least four major kinematic component 
groups over 400 km s^^ in its rest frame. HST/WFC3 images of the galaxy field show that the absorber 
is associated with two galaxies lying at 14 and 46 kpc from the QSO line of sight. The absorber is 
kinematically complex and there are no less than nine individual Mg Ilcomponents spanning 200 km 
in our Keck/HIRES optical data. The bulk of the absorbing gas traced by H I resides in two 
strong, blended component groups that possess a total logiV(H I) ~ 18 — 18.8, but most of the O VI 
is associated with two outlying components with logA^(H I) = 14.8 and 16.5. The ion ratios and 
column densities of G, N, O, Mg, Si, S, and Fe, except the O VI, can be accommodated into a simple 
photoionization model in which diffuse, low-metallicity halo gas is exposed to a photoionizing field 
from stars in the nearby galaxies that propagates into the halo at 10% efficiency. In this model, the 
clouds have neutral fractions of 1 — 10% and thus total hydrogen column densities of logAf(H) 
~ 19.5. Direct measurement of the gas metallicity is precluded by saturation of the main components 
of the H I, but we constrain the metallicity firmly within the range 0.1-1 Zq, and photoionization 
modeling indirectly indicates a subsolar metallicity of 0.05 - 0.5 Zq. This highly ionized, multiphase, 
possibly low-metallicity halo gas resembles gas with similar properties in the Milky Way halo and 
other low-redshift LLS, suggesting that at least some other galaxies have their star formation fueled 
by metal-poor gas accreting from the intergalactic medium and ionized by the stars in the host galaxy. 
As observed in the Milky Way high-velocity clouds, the strong detected O VI is not consistent with 
the photoionization scenario but is consistent with general picture in which O VI arises in interface 
material surrounding the photoionized clouds or in a hotter, diffuse component of the halo. The 
appearance of strong O VI and nine Mg II components in this system, and our review of similar 
systems in the literature, offer some support to this "interface" picture of high-velocity O VI: the 
total strength of the O VI shows a positive correlation with the number of detected components in 
the low-ionization gas. 

Subject headings: galaxies: halos, formation — quasars: absorption lines — intergalactic medium 



1. INTRODUCTION 

A full understanding of how galaxies acquire their gas 
from the IGM and return it there in the form of chemi- 
cal and kinematic feedback will likely form an important 
part of any complete picture of galaxy formation. Moti- 
vated by questions of how galaxies obtain their observed 
stellar masses and morphology, theorists have developed 
a picture in which gas enters dark matter halos and 
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galax ies by either cold flows along filaments ( Keres et al 



2005 1 , accretion of hot material that has passed through 
a shock on its way in from the IGM (Dekel & Birnboim 



two( |Maller fc Bullock| |2004). Mass loss by feedback is 
driven in "superwind" onflows fueled by m any correlated 
super novae, by radiation pressure on dust ( Murray et al. 



2006) ), or in some com plex, multiphase mixture of the 



2005 1 , or in AGN-triggered flows in those galaxies that 



possess active nuclei. Tidal or ram-pressure stripping of 
gas from satellite galaxies during gravitational encoun- 
ters may also provide a gas supply to larger galaxies, or 
at least to their halos. All these gas processes may bear 
on such observational features of galaxies as their stellar 
masses, morphologies, and colors, but to be effective they 
must act across the 100 kpc scales of galaxy halos. Unfor- 
tunately the proposed accretion, feedback, and stripping 
processes are difficult to test directly because diffuse, ion- 
ized gas in the immediate vicinity of galaxies is difficult 
to detect. We therefore have an incomplete empirical pic- 
ture of how accreting gas is distributed around galaxies of 
all types, what its temperature, density, and metallicity 
configurations look like, and how these features influence 
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Fig. 1. — The co-added COS data, plotted versus observed wavelength, and binned by 3 raw COS pixels to^ptimal sampling of 2 bins per 
resolution element. The Lyman limit system analyzed in this paper and the DLA from Mciring ct al. (2011aJ are marked at the positions 
of Lyman series lines (green), metal lines (blue), and O VI (red). The Lyman limit of the present systefnTs* clearly visible as a complete 
absorption of the QSO spectrum at ~ 1240 A. Geocoronal Lya emission has been excised near 1216 A. The solid and dashed curves mark 
the opacity of the Lyman limit system using total column densities of log A'^(H I) = 17.8 and 18.0, respectively, as described in § 3.5. 



or are influenced by the host galaxy. 

The classic quasar absorption-line technique provides 
a method for studying intergalactic and circumgalactic 
medium gas, even to very low density and metallicity. 
This technique has been used effectively to probe dif- 
fuse ionized gas in galaxy halos, over long pathlengths 
through the IGM, and throughout the halo of the Milky 
Way. The IGM samples reveal the statistical correla- 
tions of galaxies with intervening absorbers, and in some 
cases show physical properties of gas that is well within 
the virial radius. However, most such information comes 
from post-facto galaxy surveys that do not uniformly 
sample the range of galaxy properties that may be related 
to the gas properties. While extensive observational in- 
formation exists about the quantity and physical state of 
multiphase high-velocity cloud ( HVC) gas in the Milky 



Way halo ( Sembach et al. 2003; Fox et al.||2005[ [Collins 
et al. 2005'! |Shull et al.||2009| ), and gas stripping from 
dwarf galaxies is clearly demonstrated in the case of the 
Magellanic Clouds and Stream, it is not known how typ- 
ical the Milky Way is in this regard. 

To address the problem of how halo gas in other galax- 
ies is distributed and how it compares to the halo gas of 
the Milky Way, we have begun an effort to systematically 
survey the gaseous halos of low-redshift galaxies using the 
Cosmic Origins Spectrograph (COS) aboard the Hubble 
Space Telescope. Our survey chose galaxies with sight- 
lines to UV-bright QSOs passing through their halos, 
with redshifts tuned to place the AA1032,1038 doublet of 
O VI near the peak of COS sensitivity at 1250 - 1350 A. 
This survey design allows us access to all the key far-UV 
ionization diagnostics from which gas budgets, metallic- 
ities, and kinematics can be inferred. As part of this 
survey we are also obtaining galaxy spectra and high- 
resolution optical spectra of the QSOs. The main results 
of this survey will be published elsewhere; here we report 
on a serendipitous discovery from this program that nev- 



ertheless addresses its main goal: to examine the gaseous 
fuel and/or waste residing in galactic halos. 

This paper reports new HST/COS and HST/WFC3 
and Keck HIRES and LRIS data for the sightline to 
J1009-I-0713. This sightline exhibits a strong Lyman- 
limit system (LLS) bearing a wide range of ionic absorp- 
tion from multiple ionization stages and in several dis- 
tinct kinematic components. The system was discovered 
serendipitously in the data on this quasar obtained as 
part of our survey; the absorption-line system associated 
with the targeted galaxy along this sightline will be pub- 
lished separately. This paper first describes the data col- 
lection and analysis for the HST/COS and Keck/HIRES 
data on the sightline to the QSO, and on our HST/WFC3 
and Keck/LRIS data on the galaxies in the field (§ 2). 
We then report the details of our line identification and 
analysis (§p]), and on the analysis of the galaxy images 
and spectra (§[4|. Section 5 describes physical models 
for the absorber. In Section 16] we describe our general 
results and their significance ror the larger picture of gas 
in galaxy halos. We adopt the WMAP7 cosmological pa- 
rameters Ho = lOOfe = 71 km s~^ Mpc~\ = 0.734, 
and flrnh'^ = 0.1334 ( |Komatsu et"al]|201l| ). 



2. OBSERVATIONS 

2.1. COS Data 

Data on the SDSSJ100902. 06-^071343.8 sightline 
(hereafter J1009-h0713) were obtained by COS over 3 or- 
bits on 29-30 Mar 2010 as Visit 13 in Program G011598. 
Two exposures were obtained with the FUV G130M grat- 
ing at central wavelength settings 1291 and 1309, with 
total exposure times of 1497 and 2191 sec, respectively. 
Two exposures were obtained with the FUV G160M grat- 
ing at central wavelengths settings 1577 and 1600, with 
exposure times of 2002 and 2007 sec, respectively. All 
exposures were taken in TIME- TAG readout mode with 
contemporaneous wavelength stim pulses and the de- 
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Fig. 2. — The galaxy field, aligned with north at top and east to the left. The SDSS Skyserver image and our WFC3 image are shown for 
contrast. The larger panels are approximately 25 arcsec across. Galaxies 170_9 and 86_4 are at 2 = 0.356, the redshift of the LLS. At the 
redshift of the system, the 5" range bar subtends 25h~^ kpc. Precise impact parameters for the galaxies are given in Table 2. The insets 
show zoomed images of the three galaxies of interest and the QSO itself. 



fault FP-POS position. The data were processed with 
the standard CALCOS data pipeline (v2.12) on retrieval 
from the archive. More details on the performance of 
COS can be found in the GOS Instrument Handbook 
pixon et aL]|2010| and at the STScI website. 

We began our analysis with the xld files provided by 
GALGOS. These extracted ID spectra were then coadded 
to combine the multiple exposures for each grating, and 
in turn the two gratings, into a single spectrum. To align 
the four exposures as closely as possible in wavelength 
space, small (< 10 pixel) shifts were derived for each ex- 
posure individually from the profiles of strong Galactic 
ISM absorption lines. This coaddition was performed on 
the "gross counts" vector from the CALCOS xld files. 
The co-addition process sums counts from separate ex- 
posures for each pixel, keeping track of the effective expo- 
sure time for each pixel correctly. That is, if a given pixel 
is ignored in a particular exposure because of data qual- 
ity flags, that exposure time is not summed for that pixel. 
Thus we are effectively co-adding count rates rather than 
counts. Co-addition in count space allows for a simple 
derivation of the corr ect error vecto r in the Poisson limit 
of low total counts (Gehrels 1986). The final coadded 
spectrum possesses a signal-to-noise ratio S/N — 7—15 
per resolution element. The resolving power of COS is 
approximately R = A/AA = 16 - 18,000 over 1140 - 
1750 A, and the resolution element is sampled by ap- 
proximately six (analog) pixels in the ID extraction. We 
binned our final, coadded spectrum into 3-pixel bins for 
analysis. This binning leaves approximately 2 bins per 
resolution element, close to optimal sampling. 

We note that the adopted version of CALCOS did not 
apply a correction for fixed pattern noise, but that the 
coaddition of exposures taken with different central wave- 
length settings mitigates the regular shadowing pattern 
of the detector quantum-efficiency (DQE) grid wires ly- 
ing above the face of the GOS micro-channel plate de- 
tectors. During the coaddition of the data we applied a 
correction for these grid wire shadows provided to us by 



the GOS instrument team. This correction removes the 
grid-wire features, which otherwise result in ~ 15% less 
counts recorded over a few pixels and can mimic absorp- 
tion lines. 

The final, reduced COS data appear in Figure [l] The 
spectrum reveals at least two absorption line systems 
of significance along this sightline. First, we discovered 



the 



^abs 



0.356 Lyman-limit system (LLS) that is the 



main focus of this paper. This system presents an ob- 
vious Lyman-limit break at 1240 A and a host of other 
multiphase ionization stages that will be analyzed be- 
low. A damped-Lyman-a system (DLA) was discovered 
at Zgbs = 0.114 and has been analyzed by |Meiring et al.| 



( |2011a[ ) . Both of these latter systems were serendipitous 
discoveries along a sightline that was selected to pass 
through the halo of an unrelated galaxy at z < 0.25, 
and we have no reason to suspect that the selection of 
the QSO itself or the targeted galaxy introduced a bias 
in favor of these two extraordinary intervening systems. 
The properties of the targeted galaxy and its associated 
absorber will be presented as part of the main survey. 

2.2. WFC3 Data 

Based on the SDSS images, we identified two galax- 
ies that could be associated with the LLS at z = 0.356 
and/or the DLA at z = 0.114. These are labeled 
"80_5-F86_4" and "170_9" in Figure § (left panei;Q The 
indistinct profile near the QSO is partially hidden from 
view by the ground-based PSF of J1009-f0713 in the 
SDSS image. Any attempt to relate the LLS to the 
galaxy properties, or even to obtain an accurate mea- 
surement of the impact parameter of this galaxy to the 
QSO sightline, was significantly hindered by this source 
confusion. In an attempt to conclusively identify the 
LLS-associated galaxies and to possibly discover a DLA 

^ We label our galaxies with a position angle and angular sepa- 
ration in arcseconds with respect to the target QSO; thus galaxy 
170-9 has position angle of 170° degrees (N through E) and is sep- 
arated by 9" from J1009-I-0713. 
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host galaxy close to the sightline (fully within the SDSS 
PSF of J1009+0713), we obtained an image of this field 
with the Wide Field Camera 3 aboard HST. 

The WFC3/UVIS data on this field were obtained on 
24 June 2010 as Visit 44 in program G011598. We 
obtained images in two broadband filters, F390W and 
F625W. The QSO was placed at the "UVISl" aperture 
position, and the WFC3-UVIS-M0S-DITH-LINE dither 
pattern was used with an exposure time of 376 sec in 
F625W and 395 sec in F390W. This pattern resuhed in 
total exposure times of 2256 sec and 2370 sec in F625W 
and F390W, respectively. The six individual exposures in 
each filter were first processed by the default CALWFC3 
pipeline. Each exposure then had the QSO subtracted off 
using a TinyTim model of the telescope PSF before the 
six exposures were added together using MultiDrizzle. 
More details about the processing of the WFC3 im age 



can be found in the paper by jMeiring et al. (2011a) 

Even without coaddition or CjSU PSl"' subtraction, the 
WFC3 image revealed that the candidate galaxy nearest 
the QSO in the SDSS image resolves into two galaxies, 
labeled 80_5 and 86_4 in Figure J2j that are separated by 
only ~ 1". The image itself does not tell us whether 
one or more of these galaxies is associated with the DLA 
or LLS. We have attempted to constrain their redshifts 
separately with followup LRIS data, as described in the 
next section. 

2.3. Keck LRIS 

After the COS data revealed two unexpected and in- 
teresting absorbers in this sightline, we re-examined the 
imaging data from which we had selected the target QSO. 
To determine which of these galaxies, if any, were associ- 
ated with our newly detected absorbers, we obtained 1"- 
wide longslit spectra of these two galaxy candidates with 
the Keck/Low- Resolution Imaging Spectrometer (LRIS) 
on 5 April 2010. These LRIS data were taken using the 
D560 dichroic with the 600/4000 1/mm grism (blue side) 
and 600/7500 1/mm grating (red side) which gives a spec- 
tral coverage between 3000 and 5500 A (blue side), and 
5600 to 8200 A (red side). On the blue side, binning the 
data 2x2 resulted in a dispersion of 1.2 A per pixel and 
a FWHM resolution of ^-^280 km/s. On the red side, the 
data were binned 1x2, resulting in a dispersion of 2.3 
A per pixel and a FWHM resolution of ^200 km/s. Ex- 
posure times were 800s in the blue and 2 x 360 s in the 
red, which resulted in signal-to-noise ratios of at least 3 
per pixel for strong nebular emission lines in the galaxy 
spectra. 

Data reduction and calibration were carried out using 
the LowRedu?!^ IDL software package, which includes 
flat fielding to correct for pixel-to-pixel response varia- 
tions and larger scale illumination variations, wavelength 
calibration, sky subtraction, and flux calibration using 
the spectrophotometric standard star G191B2B. Precise 
and accurate systemic redshifts were obtained for both 
galaxies using a modified version of the SDSS IDL code 
"zfind" , which works by fitting smoothed template SDSS 
eigenspectra to the galaxy emission-line spectra on both 
red and blue sides. The resultant weighted-mean redshift 
for 170_9 is z = 0.355687 ± 0.00001. The long sfit used 



to obtain this first set of Keck/LRIS observations was 
oriented to run from galaxy 170_9 through the midpoint 
of the faint profile seen for 80_5-|-86_4 in the SDSS image 
(the WFC3 image was not yet available). The slit was 
1" wide, so the recorded spectrum most likely includes 
contributions from both galaxies 80_5 and 86_4. For this 
spectrum, we obtain z — 0.355574 ± 0.00002, with error 
bars calculated from statistical noise only. To account 
for systematic uncertainties in the absolute wavelength 
calibration and instrument flexure during the LRIS ex- 
posures, we adopt a larger 25 km s~^ error on the galaxy 
redshifts. The redshifts of the two galaxies are marked 
with ticks and 25 km s~^ errors at the top of Figures|3j|4j 
and[5]after translation into the rest frame of the absorber 
{zabs = 0.3558). The galaxies appear approximately 25 
km s^^ apart. Neither appears to line up exactly with 
any of the strongest absorption components, but their 
coincidence is enough to give us strong evidence of their 
kinematic association. 

Additional LRIS exposures of this field were obtained 
with LRIS in January 2011 in an attempt to separately 
constrain the redshifts of 80_5 and 86_4, which were 
blended together in the first set of LRIS exposures. LRIS 
was configured using the 400/3400 grism on the blue 
side, the 600/7500 grating on the red side, and the d650 
dichroic. This setup provided full coverage of wave- 
lengths 3000 < A < 8700 A. Three exposures, each last- 
ing 600 seconds, were obtained in ^ 1 arcsecond seeing. 
Each exposure had a different position angle, with the 
goal of obtaining spectra of 80_5, 86_4, and other sources, 
while minimizing contaminating flux from the QSO. This 
observation confirmed the detection of the [O III] and II/3 
lines at z = 0.355 for 86_4. With the separation of 80_5 
and 86_4 measured from the WFC3 image, we were able 
to marginally separate the spectral traces of 80_5 and 
86_4. We do not find the same pattern of emission lines 
at the position of 80_5, so we conclude that this galaxy 
is not at the same redshift as 86_4. The best solution for 
the weak emission line detections at this position gives a 
solution of z — 0.879. We therefore disregard this galaxy 
in our subsequent analysis. 

2.4. Keck HIRES 

Using the Keck I High Resolution Echelle Spectrom- 
eter with the blue collimator (HIRESb), we obtained a 
high dispersion spectrum of the J1009-I-0713 on 26 March 
2010. These data were taken through the CI decker, re- 
sulting in a FWHM resolution of « 6 km s^^. The grat- 
ing angles were set to provide wavelength coverage from 
3050 - 5870 A, with two small gaps related to the sepa- 
ration of the three CCD mosaic. We integrated for two 
1800s exposures under good conditions. The data were 
reduced with the HIRedux pipeline]^ to flatfield, sky sub- 
tract, wavelength calibrate, and extract the spectra. The 
data were optimally coadded and normalized to unit flux 
by fitting a series of Chebyshev polynomials. The final 
spectrum has a S/N of 15 per 1.3 km s~^ pixel redwards 
of 3800A decreasing to S/N=5 at 3200 A. 

3. ABSORPTION LINE IDENTIFICATION AND ANALYSIS 
3.1. General Approach 



® http:/ /www. ucolick.org/~xavier/LowRedux/index. html 



^ http:/ /www. ucolick.org/~xavier/HIRedux/index. html 
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Velocity [km s '] 

Fig. 3. — Lyman series lines for the z = 0.355 LLS. The velocity 
zero point is set to z = 0.3558. The points with error bars at the 
top mark the measured velocities of the galaxies with 25 km s~^ 
uncertainty. 

To identify and fit absorption lines in the spectrum of 
J1009+0713 we have adopted the vacu um wave l ength s 
and atomic data from the compilation by Morton ( 2003 ) . 
Line measurements appear in Table [l] Most column den- 
sities were obtained by direct integration over the line 
profile, except in the case of the optical lines, some satu- 
rated UV lines, and the Lyman series lines in the outlying 
components, where profile fits were used. 

The line profile fits use a custom Voigt-profile fitting 
software that attempts to optimize a model of the line 
column densities, doppler b parameters, and velocities for 
an arbitrary input set of lines and atomic data, by min- 
imizing the x^- This cod e was first used to analy ze the 
PG1211-M43 sightline by|Tumlinson et al.|(|2005|. This 



software first constructs intrinsic line profiles and then 
convolves these with the appropriate non-Gaussian COS 



line-spread f unction (LSF) as described by Ghavamian 
et al.| ( |2009[ ) and available on the STScI COS website 
These tabulated LSFs are specified for 50 A intervals 
through the range of the G130M and G160M gratings. 
For line fits we adopt the nearest LSF grid point based 
on the observed wavelength of the lines being fitted. 

3.2. The Lyman Limit System 

The z = 0.3558 absorber exhibits a line-black Lyman 
limit at 1236 A in the observed frame. This system also 
appears in at least 14 distinct Lyman series transitions, 
the first twelve of which appear in Figure[3j Inspection of 
the strong Lya profile reveals absorption extending over 
^ 400 km s~^. No damping wings are apparent, which 
limits the total column density to logiV(H I) < 19. At 
Ly/3 - e, a separate component at -1-180 km s""'^ becomes 
distinct, which we term component D. By Lyry, a sep- 
arate component at —95 km s^^ is visible and distinct 
until blending between the Lyman series lines themselves 
becomes severe at Ly/.t. The two innermost components 
labeled B and C in Figure [3] never separate in the Ly- 
man series lines and are defined by their distinct profiles 
in the metal lines, particularly Mg 11, as shown below. 
Component groups A, B, and C all break further into 
multiple blended components in the HIRES data. 

We note that these components are defined, wherever 
possible, by the H I absorption or by the metal lines, as 
noted in Table 1. The number of apparent components 
is larger in the higher resolution optical data covering 
Mg II and Fe II. This finding suggests that there may be 
further component structure in the H I and metal lines 
that we cannot resolve at the R — 18,000 resolution of 
COS. In what follows, we analyze the major component 
groups as coherent objects with the proviso that they 
may in fact consistent of subcomponents below the level 
of our spectral resolution. 

3.3. Metal lines and ionization 

The COS data on the LLS reveal absorption from a 
wide range of ionization stages, as shown in Figure |4j 
We detect absorption by C, N, O, Mg, Ca, Si, S, and Fe 
in the COS and/or HIRES data, and in ionization stages 
from Mg I to O VI. However, the overall impression is of 
an absorption-line system that possesses a high degree 
of ionization. The species Mg I and O I are the only 
neutral atoms that are convincingly detected; Ca II is 
usually detected in mostly neutral gas but is very weak 
here. Limits are placed on the neutrals of C, N, Si, and 
Fe. By contrast, strong absorption is detected in C III, 
N III, Si III, and Fe III. Our COS data do not cover C IV 
or Si IV, but N V AA1238, 1242, typically weak, is not 
detected in the S /N ^ 5 data at that wavelength. These 
detections and limits, apart from any detailed analysis, 
point to gas that is substantially ionized in all compo- 
nents. 

3.4. The VI Absorption 

The strength and kinematics of the O VI in this ab- 
sorber are extraordinary. Integrating over the full pro- 
file, this is among the strongest intervening intergalactic 
O VI absorbers yet detected, with a rest-frame equiva- 
lent width Wr = 835 ± 49 mA in the A1032 profile (cf. 
Tripp et al.|2008| [Thorn fc Chen|[2008l [Danforth fc Shull 
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Fig. 4. — Metal absorption lines from the z = 0.3558 LLS. Vertical dashed lines mark the centroid velocities of components A - D. 
Transitions are ordered by atomic element and, within that, by increasing ionization potential. The points with error bars at the top mark 
the measured velocities of the two associated galaxies with 25 km s~^ uncertainty. 
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TABLE 1 

Measured Component Column Densities in the z = 0.3558 System 



Line 


A 


B 


C 


D 


Centroid 


— 95 km s""*^ 


— 10 km 


+60 km s^^ 


+180 km 


Range 


-200 to -65 


-65 to 30 


30 to 130 


130 to 250 


Lya 1215 


blend 


blend 


blend 


Ly/3 1025 


blend 


blend 


14.83,+183,17« 


T ■\7>~^J Q79 


blend 


blend 


1 A 7Q -ui an 1 fi^ 


Ly5 949 


blend 


blend 


contaminated 


Lye 937 


blend 


blend 


contaminated 


LyC 930 


blend 


blend 




Ijvn 926 
y I 




blend 




Lye 923 


16.33,-93,22=^ 


blend 




Lyt 920 


contaminated 


blend 




LyK 919 


16.44,-94,20^ 


blend 




LyA 918 


not fitted^ 


blend 




Ly^ 917 


blend 


blend 




C I 1157 


< 14.4 


< 14.3 


< 14.3 


< 14.4 


C II 1036'^ 


14.0 ±0.1 


14.6 ± 0.3 


14.6 ± 0.3 


< 13.5 


C III 977^ 


> 14.3 ±0.2 


> 14.3 ±0.2 


> 14.2 ±0.2 


> 13.5 + 0.1 


N I 1199, 1200 


< 13.7 


< 13.7 


< 13.7 


< 13.7 


N II 1083 


< 13.4 


14.1 ± 0.1 


14.0 ±0.1 


< 13.6 


N III 989'' 


14.3 ±0.1 


14.6 ±0.1 


14.5 ±0.1 


14.5 ±0.1 


1S.T A/ 1 OQC 

IN V izoo, iZ4tZ 




< 13.7 


< 13.7 


•C lo. ( 


O I 1302 


< 13.9 


14.6 ± 0.2 


14.1 ± 0.2 


< 13.9 


O I 988 


< 14.0 


14.4 ± 0.2 


14.2 ± 0.2 


C 


OVI 1032 


14.67 ±0.1 


14.40 ±0.1 


14.40 ±0.1 


14.29 ±0.1 


OVI 1038 


14.56 ± 0.1 


1 4 QQ -L n 1 


14 26 + n 1 


14.33 ± 0.1 


Mg I 2852 


< 11.2 


11.9 ±0.04 


11.4 ±0.1 


< 11.2 


TV J TT o ^ R 

Mg 11 2796 


12.82 ± 0.02 


13.59 ±0.04 


13.45 ± 0.04 


< 11.2 


Mg il zoUo 


1 o no _L rx no 

12. yy ± U.U2 


13.75 ± 0.03 


13.62 ± 0.03 


< 11.5 


Si I 2515 


< 11.7 


< 11.7 


< 11.7 


< 11.7 


Si II 1260* 


> 12.7 lb 0.3 


> 13.6 ±0.2 


> 13.6 ±0.2 


< 12.7 


Si II 1190 


c 


13.9 ±0.2 


13.8 ±0.2 


< 13.6 


Si II 1193 


... d 


13.8 ±0.2 


13.5 ±0.2 


< 13.4 


Qi TTT 1 Oflfif 
Ol 111 iZUO 




> 13.8 ±0.2 


> 13.8 ±0.2 


\ 1 o -U n o 
>■ iz.o it yj.z 


S II 1253S 


< 14.8 


< 14.7 


< 14.7 


< 14.7 


S III 1012 


14.2 ±0.1 


14.3 ±0.1 


13.8 ±0.1 


< 13.4 


S IV 1062 


14.0 ±0.1 


< 13.8 


< 13.8 


< 13.8 


Ca II 3934 


< 11.4 


11.50 ±0.15 


< 11.4 


< 11.4 


Ca II 3968 


< 11.9 


< 11.9 


< 11.9 


< 11.9 


T^i TT QQO/I 

1 1 ii ooo4 


< 11.5 


< 11.5 


< 11.5 


< 11.5 


Fe I 2484 


< 11.3 


< 11.3 


< 11.3 


< 11.3 


Fe II 2382 


< 11.8 


13.71 ±0.4 


13.45 ±0.16 


< 11.8 


Fe II 2600 


< 11.6 


13.65 ± 0.07 


13.37 ± 0.02 


< 11.6 


Fe II 2586 


< 12.2 


13.88 ± 0.03 


13.38 ± 0.04 


< 12.2 


Fe II 1144 


< 13.8 


14.1 ±0.2 


... c 


... c 


Fe II 2374 


< 12.8 


14.02 ± 0.09 


< 12.8 


< 12.8 


Fe II 1063 


< 13.8 


13.7 ±0.2 


13.6 ±0.2 


< 13.8 


Fe III 1122 


... c 


14.3 ±0.1 


14.2 ±0.1 


< 13.9 



Note. — All column density measurements were obtained via direct integra- 
tion of the line profiles over the given velocity ranges, except where noted in the 
comments field. Upper limits for UV lines are obtained as 2a limits by direct inte- 
gration of the apparent column densities over the stated component velocity ranges. 
For multiplets, the limits use the strongest uncontaminated transition. Line pro- 
file fitting results are given as a triplet; column density, velocity centroid, doppler 
b- value. 

^ Part of the definition of this component; HI where possible, metal lines where 
necessary. 

^ N III A989.80 column densities are uncertain because they are blended with Si II 
A989.87 to an unknown degree. These column densities should be treated as upper 
limits. Component D is even more uncertain, because it appears to be blended 
with other unidentified absorption. 

Contaminated by unknown absorption - no measurement or limit is possible. 

Very noisy. 
" Subject to mild saturation. 

' Strongly saturated and /or line black absorption. The equivalent widths over each 
velocity range are converted to column densities assuming a linear curve of growth, 
which yield robust lower limits. 

s S II A1259 is contaminated and A1253 and A1250 are undetected. The A1253 line 
provides the best upper limits. 
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Fig. 5. — Kinematic structure in apparent column densities ( |Sav-| 
|age fc Sembach|1991| l for low ions and O VI. Major components A, 
a, and (J are broken down further according to their appearance 
in the Mg II and Fe II data obtained by HIRES. Note the wide 
range of apparent column density ratios between the low ions and 
O VI, particularly in components A and D. Flat tops arise from 
the truncation of the flux at 0.01 in the normalized data to avoid 
unphysical values in the optical depth for saturated lines. 



2008[|^° As shown in Figure [4j this absorption spreads 
over the fuU 400 km s^^ range of the system, in a nearly 
flat, complex profile. Figure^ shows the O VI profile af- 
ter it has been conv erted into apparent colum n densities 
using the method of Savage & Sembach ( 1991 ). Account- 
ing for a modest degree of noise in the data, the A1032 
and A1038 profiles agree well over the full velocity range, 
supporting the case that all this absorption is truly O VI, 
and that it covers the full kinematic range of absorption 
by the other detected ionization stages. While compo- 
nent groups A and D appear distinctly in the O VI pro- 
files, B and C are not clearly separated; they are both 
blended together and possibly too broad to show indi- 
vidual peaks in the column density plot (Figure [5]). We 
conclude that O VI is associated with all four identified 
component groups in this system. 

Another notable feature of the O VI is its relationship 
to the neutral and low-ionization gas in the various com- 
ponents. In component A, strong O VI at logiV(0 VI) 

The reported values for some of the other Wr > 300 mA 
systems vary from study to study, depending on how each system 
was broken down into components. Our reported value takes in all 
the identifiable components. This pattern emphasizes the difficulty 
of interpreting kinematically complex absorbers in terms of simple 
models, which we will take up below. 



= 14.6 coincides with strong C III, N III, and Si III, but 
also Mg II and C II. In components B and C, somewhat 
weaker O VI is associated with much stronger low ions, 
such as C II, Mg II, Si II, and Fe II. In component D, 
only H I, C III, N III, Si III, and O VI are detected, 
with no corresponding Mg II. Clearly the state of ioniza- 
tion varies across the absorber, which has an intrinsically 
multiphase structure. 

3.5. HI Solution 

With the four components defined either by their H I 
or metal- line profiles (see notes to Table jTl), we can now 
derive a solution for the H I column densities in this 
absorber. Components A and D both appear in 3 or 4 
higher Lyman series profiles, and can be profile-fitted di- 
rectly. However, both B and C are strongly saturated 
and blended together, so simple line-profile fitting can- 
not disentangle them and measure their column densi- 
ties without additional assumptions. The two straight- 
forward components are taken up first. All relative ve- 
locities are given with respect to the adopted systemic 
redshift of Zabs = 0.3558. 

Component A: This component is defined at —95 km 
s~^ by Lyyy, 6*, k, and A, where it appears distinctly but 
is blended with the blue wing of component B. It ap- 
pears to be contaminated with unrelated absorption at 
Lyi. Profile fits to Lyr;, 9, and k were attempted, with 
the results appearing in Table [l] Consistent results are 
obtained with log7V(H I) = 16.42 ± 0.08, 6 = 20 ± 3 km 
s~^, at —94 ± 1 km s^^. These values are adopted for 
further analysis. 

Component D: This component is defined at -1-180 km 
s~^ by Ly/3, 7, and perhaps 5. It is inseparably blended 
with component C of Lya, contaminated at Lye, and 
absent in higher Lyman series lines. Only Ly/3 and Ly7 
give acceptable fits, with log N(II I) = 14.8, 6 = 16 ± 3 
km s^^ at 180 km s^^. This fit is adopted for further 
analysis. 

Components B and C: These two components are both 
strongly saturated and severely blended with one another 
at a velocity separation of 60-80 km s~^ (as shown by 
their metal lines) . Because of this blending and satura- 
tion we cannot, even in principle, fit them and constrain 
their column densities separately. Instead, we attempt 
to impose firm limits on the total NiYi I), based on the 
Lyman limit opacity and observed Lyman series line pro- 
files. 

A firm lower limit to the total iV(II I) comes from the 
observed Lyman-limit opacity. For closely-spaced com- 
ponents, this limit is insensitive to the separations and 
column density ratios and traces mainly the total column 
density of H I. We display two limiting cases in Figure [T 
We adopt the Lyma n continuum opacit y as function of 
wavelength given by [Osterbro"^ ( |1989[ ). The red solid 
curve shows log7V(H I) = 17.8, which leaves ^ 5% trans- 
mission across the detected range of the Lyman limit sys- 
tem. By log7V(II I) — 18.0, the transmission is no longer 
visible and matches the absence of source flux in this re- 
gion of the spectrum; higher column densities are also 
consistent with the data. We therefore regard log 7V(H I) 
= 17.8 as a firm, conservative lower limit to the total col- 
umn density of the LLS. Given the Lyman-limit opacity, 
17.8 seems to have a low probability, so we also adopt 
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Fig. 6. — The LLS Lya profile with model lines overlaid, showing 
the effects of line saturation. The top panel shows a model with 
logAf(H I) = 18.5 in each of the components B and C, for a total 
of log A'^(H I) = 18.8. The combined profile, shown with the heavy 
black line, is still a marginally acceptable fit to the data. The 
lower panel increases the individual components to logAf(H I) = 
18.7 and the total to 19.0; by this point damping wings no longer 
fit the data. 

a slightly higher value of log iV(H I) to bound the lower 
end of a 95% confidence interval. 

Upper limits to total 7V(H I) come from the absence 
of damping wings from the Lya profile of the LLS. This 
limit is somewhat less firm than the limit from continuum 
opacity, since it varies slightly with the line broadening 
and number of components. Two identical components 
at the positions of the B and C components with b = 20 
km s~^ and log A^(H I) = 18.5 show damping wings that 
are only marginally consistent with the data (Figure [6]). 
When each component is increased to logiV(H I) = 18T7 
(for a total of 19.0), the damping wings are clearly ex- 
cluded by the data. These profiles and limits strictly 
depend on the number of assumed components, which 
we can infer from the HIRES data. If we assume com- 
ponents at each of the three strongest Mg I and Mg II 
components observed in the HIRES optical data, which 
are likely to trace the majority of the total H I, we 
find that by a total of logA^(H I) = 18.5, the damping 
wings appear, and higher values are increasingly poor 
fits. This is consistent with the two-component analysis. 
As the number of assumed components within the inner 
80 km increases, the column density per component 
increases. By log A^(H I) = 19, no tuning of b and number 
of components (even one) can hide the damping wings, 
so we regard this value as a firm upper limit to the total 
N(H I). The weak presence of damping wings at a total 
logAr(H I) — 18.8 suggest that the value is likely to be 
lower than that, so we adopt this as an upper bound to a 
95% confidence interval. For further analysis, we adopt 
at 95% confidence interval of logA^(H I) = 18.0 - 18.8, 
with a flat probability density within that range, since 
we have no information about which values in this range 
are more likely. 

4. METALLICITY AND STAR FORMATION OF THE 
NEARBY GALAXIES 

The Keck/LRIS and WFC3/UVIS data on the galax- 
ies in the field provide useful information about their 
redshift, luminosity, star formation rate, and metallic- 
ity. The reduced Keck/LRIS spectra appear in Figure [t] 
and the calibrated and PSF-subtracted images appear in 
Figure [2j 



At z = 0.3558, [OH] AA 3727, H/3, and 
[OIII]A4959A5007 are readily detected. The strengths of 
these emission lines, with appropriate corrections, pro- 
vide constraints on th e star formation rate and metal- 
licity of the galaxies ( [Levesque et al. 2010). Accord- 
ingly, the accuracy of these quantities depends upon the 
accuracy of the relative emission line strengths and on 
the degree to which the flux-calibration of the spectra is 
"absolute." We therefore try to ensure that the flux cal- 
ibrations of the blue and red sides of the LRIS spectra 
are consistent with each other, and attempt to correct 
for light losses in the 1" slit by comparing our spectra 
with SDSS photometry if it is available for our galax- 
ies. We accomplish this absolute flux correction by con- 
volving the LRIS spe ctra with SDSS ugriz filters (see 



Da Silva et al. (2011) for a description of the IDL code 
"spec2mag") arid comparing the observed spectral ap- 
parent magnitudes with the SDSS catalog apparent mag- 
nitudes (see Werk et al., in prep, for details). 

Only galaxy 170_9 has available SDSS photometry that 
enables us to perform an absolute flux correction. We do 
not correct the spectral fluxes of objects 80_5 and 86_4 
since they do not have SDSS photometry (they are not 
even detected as one galaxy). The resultant corrective 
flux factor for 170_9 is 1.58 (0.494 magnitudes in the 
g— band) on the blue side and 1.65 (0.543 magnitudes 
in the i— band) on the red side. Since we have ugriz 
photometry available to us for ob ject 170_9, we are able 
to use the IDL-code kcorrect (Blanton & Roweis 2007 



v4_2) to obtain an estimate of its stellar mass, approxi 
mately 5.5x10^ M©. 

We estimate the SF R from the H/3 emissio n line flux 
using the calibration of Calzetti et al. (2010), such that 
SFR (H;3) = SFR (Ha) / '2M, where^FR (Ha) (Mq 
yr-i) = 5.45 X 10-42 ^^.^ g-i^ ^he factor of 2.86 

represents the intrinsic ratio of Ha to H/3 for case B re- 
combination at an effective ter nperature of 10,000 K and 
electron density of 100 cm"'' (Hummer & Storey| 
The SFRs appear in Table 2. 



1987). 

i-'or 86_4, the Sl<'k is a 



formal lower limit, since we assume there was significant 
light loss from using a 1" slit (on the order of ~38%). 
The errors given in the table account for only the RMS of 
the line measurements; the flux calibration uncertainty, 
read noise, sky noise, and flat-fielding errors, are approx- 
imately ^-^1% for the emission lines in 170_9 and ~7% for 
the lower S/N emission lines in 86_4. 

We obtain the nebular oxygen abundances for galaxies 
86_4 and 170_9 usi ng the stron g line R23 method origi- 



nally presented by Pagel et al. 



calibration of McGaugh ( 1991 ) 



[([OH] AA3727 Tpi 



(1979), according to the 
. R23 is defined as log 
OIII] A5007)/H/3]. The 



drawbacks of the R23 method include a well known de- 
generacy and turnover at ^ O.3Z0 and large systematic 
errors on the order of 0.25 dex due to ag e effects and 



stellar distributions (Ercolano et al. 2007). To correct 



the emission lines for foreground reddening, we use a 
reddening functio n normalized a t H/3 f rom the Galactic 
reddening law of 



A„/E(B-V) = 3. 



C ardelli et al 
and using an 



(1989) assuming R„ 

^ „ E(B=V) value of 0.013 

( Schlegel et al.|1998 |. We cannot break the R23 degener- 
acy using any of the known methods (e.g. [Nil] /[OH]) be- 
cause our LRIS data do not cover lines redward of ^6500 
A in the rest frame of the galaxies. The upper and lower 
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Rest-Frame Wavelength [ A ] 
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TABLE 2 
Galaxy Properties 



] 4900 5000 5100 5200 5300 6400 6500 6600 6700 6800 6900 

Observed Wavelength [ A ] 
Rest-Frame Wavelength [ A ] 

3600 3700 3800 3900 4800 4900 5000 5100 



[Oil] 



Hp [Ollll 



4300 4900 5000 5100 5200 5300 6400 6500 6600 6700 6800 6900 

Observed Wavelength [ A ] 

Fig. 7. — Keck/LRIS spectra of galaxies in the field for 86_4 
(top), and 170-9 (bottom). The break in the wavelength range is 
placed to avoid LRIS dichroic near 5500 A. Detected emission lines 
are labeled. The red and blue sides for 86_4 are not normalized as 
they are for 170_9 (see § 2.3). 

branch metallicities are presented in Table 2. Based 
on the well-known ma ss-metallicity relation of galaxies 
(Tremonti et al. 2004), it is reasonable to assume the 
upper- branch metallicity for 170_9 given its estimated 
stellar mass of 5.5 x 10^ M©. The upper branch value 
is consistent with solar metallicity according to the up- 
dated solar oxygen abundance of 12 -|- log(0/g) = 8.69 



(Allende Prieto et al. 2001 Asplund et al. 2009). The 



spectrum ot 8b_4 mdicates a slightly lower upper-branch 
metallicity than the brighter galaxy 170_9, but given the 
errors in the measurements, all the upper branch mea- 
surements are consistent with solar metallicity and with 
one another. 

As described below, we unfortunately have no di- 
rect measurement of the gas metallicity in the strongest 
components of the LLS. However, the indirectly indi- 
cated metallicity from photoionization modeling is Z = 
0.05 — 0.5^0. If this is correct, the gas would appear 
to have a low metallicity with respect to the galaxies 
and would favor an explanation in which the gas has en- 
tered this system from some other source - the IGM or 
a stripped dwarf galaxy. However, given the uncertainty 
in the metallicity of the strong LLS components, and 
the upper/lower branch ambiguity for the galaxies, we 
cannot draw firm conclusions from this metallicity com- 
parison. 



Property 


80_5 


86_4 


170_9 


Rcdshift'^ 


0.88(1) 


0.355574(2) 


0.355687(1 


P [kpc] 




14.25 


46.48 




23.1 ± 0.2 


24.1 ± 0.2 


22.2 ±0.1 


"16251^ 


22.5 ±0.1 


23.2 ±0.2 


21.3 ±0.1 


[O/H] (upper)<^ 




-0.2 


0.0 


[0/H] (lower) 




-0.6 


-0.9 


SFR [Mq yr-l] 




> 0.2 


2.1 ± 0.02 



To account for systematic effects in addition to formal 
fitting errors, we adopt a ±25 km s~^ uncertainty for these 
redshifts. 

^ Broadb and AB magni t udes d erived from the imaging 
analysis of |Meiring et al!] l |2011a| . 

These estimates are equivaleiit to [O/H] = log (O/H) 
— log (0/H)q, where the solar abundance corres ponds to 
[O/H] = or 12 + lo gfO/H) = 8.69 (Allende Prie to et al.| 
|2001| Asplund et al.|2009| ). These estimates carry an error 
oi ±6 lb dex. 



5. PHYSICAL MODELING AND INTERPRETATION 

In this section we move from measuring column den- 
sities and kinematics to interpreting these quantities in 
terms of physical models, focusing on what we can learn 
about the ionization, metallicity, and total gas budgets 
in these components from the detected absorption. 

From the detections of the various ionization stages of 
C, N, O, Mg, S, Si, and Fe alone it is evident that at least 
some of the gas in this absorber possesses a high degree of 
ionization, even in the higher column density component 
groups B and C. These two components exhibit absorp- 
tion in the second ions of C, N, S, Si, and Fe, the third 
ion of S, O VI, and relatively little absorption from the 
neutrals and first ions of these elements. Apart from any 
detailed modeling of ionization, the ratios N II / N 1 and 
Fe 111 / Fe II indicate a significant ionization fraction. 
The solid detection of O VI in all components further 
indicates a highly ionized fraction of the gas, whether or 
or not it arises from the same ionization mechanism as 
the lower ions. 

To estimate the degree of ionization in the absorb- 
ing clouds and from that the total budget of gas that 
has been detected, we produced photoionzation mod- 
els with the Cloudy modeling cod e (version 08.00) last 
described by Ferland et al. (19981. This effort is hin- 
dered by several limitations that must be acknowledged 
at the outset. Photoionization models described by a 
density, ionization parameter, incident spectrum, metal- 
licity, and total gas content require a larger number of 
uncertain parameters than are formally constrained by 
the data, so any results must be considered indicative 
rather than conclusive. The shape and orientation of the 
absorbing cloud with respect to the nearby galaxies and 
to the observer are both unknown; we can only assume 
the simplest cases. While these two factors are generic 
to absorption systems in QSO sightlines, this particular 
system presents two additional features that complicate 
modeling. Its two strongest components, labeled B and 
C, are inseparably blended with one another in all avail- 
able lines of H 1, so models are not robustly constrained 
as usual by a definite N(H I); we have only a loose con- 
straint on the sum of the neutral hydrogen column in 
B and C (see § 3.5). Furthermore, many of the most 
important lines of intermediate ions, such as C III and 
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Si III, are strongly saturated and blended across all the 
component groups. 

With these caveats in mind, we model the absorb- 
ing clouds with a single plane-parallel slab model illu- 
minated on one side by a uniform incident radiation 
field. The model field uses two components. First, we 
attempted models with the composite extragalactic radi- 
ation field of Haardt & Madau (HMOS, as implemented 
within Cloudy). This spectrum yields ^hm = 3.2 x 10** 
photons cm for > 1 Ryd photons at each point in 
space. By itself, the composite HM field generally fails 
to achieve the observed ratios of second to first ions of 
C, N, O, Si, and Fe until the ionization parameter is in- 
creased by a low density and the clouds are more than 
1 Mpc in size, too large to reside in galaxy halos. It is 
perhaps not surprising that this spectrum fails, since it 
represents the average ionizing background of all galaxies 
and QSOs as transmitted by the IGM. But our absorbers 
are not located at a random point in space, but within 15 
- 50 kpc of galaxies with star formation rates of ~ 0.2 — 2 



M, 



yr 



Models of the Milky Way halo that attempt to 



account for the ionization of the Galactic high-velocity 
clouds typically find that the Milky Way may contribute 
10^ ionizing photons cm~^ s~^ even ~ 100 kpc into the 



halo in the direction normal to the d isk (Bland-Hawthorn 
fc Maloney|l999b|a[[fb^ e t al.|2005[), and nearly 10" pho- 
tons cm~^ s~^ at 40 kpcf^ This flux exceeds the con- 
tribution of ^HM by a significant margin. Thus even 
a modest degree of star formation, if its UV radiation 
propagates efficiently into the halo, can exceed the ex- 
tragalactic background at physical separations like those 
observed in the present system. 

We therefore add in a second ionizing field compo- 
nent, modeled as a galaxy undergoing a recent episode 
of star formation, drawn from the Starburst99 libraries 



(Leitherer et al. 1999) and corresponding to a continu- 
ous star-formation history over 10 Myr with a Salpeter 
IMF from 1 - 100 AIq and solar metallicity (the ionizing 
photon output at 0.5^0 is only 15% higher). For a star 
formation history extending far into the past, the instan- 
taneous spectrum reaches a steady state beyond a few 
million years and changes negligibly with the assumed 
age. The galaxy metallicities (see § 4) suggest that we 
should adopt Zq for the starburst. This spectrum gen- 
erates a total of = lO^S " (SFR / Mq yr'^) ionizing 
photons s^^. We assume that a uniform, wavelength- 
independent fraction f^sc of these ionizing photons es- 
cape the ISM of the galaxy to illuminate material in the 
halo. At a distance d in kpc, this ionizing field will impart 
to a cloud of uniform density ri/f a field characterized by 
the dimensionless ionization parameter U: 



U = 



1 



Q*fesc 



HM 



0.01 cm-3 



10 



-3 SFR fesc 

Mqyt-^ oT 



20 kpc 
d 



10" 



where in the final equation the quantities have been 

These models arc consiste nt with the Ha emission from the 
HVCs fWeiner et al. 2002 Putma n et al.|2003[ l, but not the Mag- 
ellanic Stream, which is hypothesized to have at least some of its 
gas ionized by shocks ^Bland-Hawthorn et al._2007] l. 



scaled to typical values for halo clouds. This relation 
does not account for radiative transfer in the Lyman con- 
tinuum, and assumes normally incident radiation, and so 
it is only an approximation for real clouds. For d = 20 
kpc and SFR — 1 Mq yr~^, the ionization parameter 
\ogU ^ —3. The constants 10"'^ and 10"'' are conve- 
nient approximations good to better than 10%. In these 
typical conditions, the stellar ionizing fiux exceeds the 
HM spectrum by a order of magnitude, consistent with 
results from the Milky Way. At all higher SFRs, higher 
fescj a-nd lower distances, the stellar ionizing field is even 
more dominant. Since we do not know the cloud dis- 
tance to the galaxy, the density, or the fesc, we cannot 
use this relation to specify U for a given absorber. How- 
ever, once a photoionization model has been found and 
U constrained by observations, along with SFR for the 
galaxies, we can use this relation to obtain estimates of 
the other parameters. 

The sum of these HM05 and starburst radiation field 
components is normalized within Cloudy to achieve a 
specific ionization parameter at the face of the model 
cloud. The model clouds are calculated from their 
sharp illuminated face to a stopping column density of 
logA^(H I) = 19.5, to cover components B and C; these 
models are valid for any lower total A^(H I) if their results 
are truncated at that column density so are applicable to 
A and D as well. The model clouds are assumed to have 
a uniform metallicity of appropriately scaled solar rela- 
tive abundances; models with total metallicity from solar 
to 1% solar were explored. 

5.1. Components B and C 
5.1.1. Metallicity 

In B and C we detect weak absorption from neutral 
oxygen in two transitions: O I A1302 and A988. These 
lines are likely unsaturated, as the apparent column den- 
sity profiles shown in Figure [5] (middle panel) are con- 
sistent in the strong component groups B and C. The 
B component of the A 1302 line appears mildly contami- 
nated; we adopt the value logA^(0 I)= 14.5 ±0.2, boost- 
ing the error from 0.1 to 0.2 dex to account for the uncer- 
tainty in the contamination. For component C, we ob- 
tain consistent measurements from both lines, log N{0 I) 
= 14.2 ± 0.2. Significant absorption is not detected in 
components A or D, and direct integrations over their 
velocity ranges yield upper limits of logAf(OI) < 13.9 
(2^)- 

Neutral oxygen is a useful metallicity indicator in mod- 
erately ionized gas, since its ionization potential of 13.6 
eV and charge exchange reactions with H lock its ioniza- 
tion fraction to that of hydrogen in typical conditions. 
Unfortunately these well-detected O I lines must be com- 
pared to a highly uncertain total H I column, ranging 
anywhere from log iV(H I) = 18.0 - 18.8 (95% confidence). 
With this range, metallicities from 0.1 — I.QZq are all 
permitted in comparison with the updated solar oxygen 
abundance 12 + log(0/g ) = 8.69 (|Allende Prieto et al. 



2001 Asplund et al.|2009 ). We therefore cannot obtain a 
robust, direct measurement of gas metallicity in this sys- 
tem, owing to the systematic error in A^(H I) introduced 
by the strong saturation and blending of the two inner 
components. However, we can impose upper and lower 
limits, Z = 0.15 — I.QZq based on the 95% confidence 
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Fig. 8. — Photoionization modeling for the B and C velocity com- 
ponents of the LLS. This model assumes the Starburst99 spectrum 
described in the text with log (7 = —3.5 and Z = O.IZq. The left 
column shows column density ratios of key ions, while the right 
shows the cumulative column densities. Both are plotted with re- 
spect to cumulative N(H I). The ranges permitted by the data for 
column densities and ratios are marked, coded in color to match 
their model curve. In the right column, the species are color coded 
according to their stage of ionization, e.g. 1 equals C I in orange, 
2 equals C II in blue. 

range derived for A^(H I) above (or 0.1 — I.GZq with the 
firm hmits on iV(H I)). This circumstance somewhat hm- 
its our abihty to interpret this material as either infall 
or outflow from the associated galaxies. However, mod- 
eling of absorption from ionized gas indirectly indicates 
sub-solar metallicity, as shown in the next section. 

5.1.2. Ionization 

The results of our Cloudy photoionization model with 
this Starburst99 input spectrum appear in Figure |8l 
specifically for the parameters \ogU — —3.5, and 10% 
solar metallicity in the gas. The left column of this fig- 
ure shows the column density ratios of the key ions, and 
the right column shows the cumulative column densi- 
ties, iV(X), as a function of H I column density into 
the cloud (note that this is not physical depth, or total 
H column, which is proportional to depth for constant 
density). These quantities are plotted with respect to a 
varying 7V(H 1) since we must account for the uncertain 



total H I column density. If that column density were 
well-constrained by measurement, we could choose one 
vertical slice though this space with a fixed A^(H I) and 
analyze models of varying U and Z, as is typically done 
for photoionization models. Instead, we analyze these 
families of models with a range of A^(H I) and proceed 
in two steps. First, we adjust the ionization parame- 
ter U to place the ion ratios in the left column within 
their permitted ranges. Then, holding U fixed, we vary 
the metallicity until the cumulative column densities in 
the right column fall into the observed ranges for them- 
selves individually and also for H I. Where these curves 
lie inside the boxes that mark data, both the total N(X) 
and the A^(H 1) are reproduced by that particular model. 
Thus, the ion ratios constrain the ionization parameter 
while the column densities constrain the metallicity. 

The model with log J7 — —3.5 provides an acceptable 
fit to the ratios of the observed ions, except that it un- 
derproduces N III and produces essentially no O VI or 
S IV. However, the detected N III is likely contaminated 
by Si II or other absorption and so the observed ratio 
N III / N II is artificially high; we discount this ratio to 
discriminate models. Otherwise, these parameters are a 
good model for the components B and C. At the depth 
into the cloud at which this model has achieved A^(H I) 
~ 18 — 18.8 the clouds are still predominantly ionized. 
As shown in Figure 8, the lower A^(H I) models within 
the permitted range have lower neutral fractions, while 
the higher iV(H I) have higher neutral fractions, so that 
the typical total N(H) range of these models is Af(H) 
~ 19.5 — 20. Thus in models with pure photoionization, 
these clouds are nearly all ionized. 

The column densities of well-measured ionization 
stages (right column) also indicate a subsolar metallic- 
ity. This constraint on the metallicity should be consid- 
ered indirect, because it depends on the assumption of 
pure photoionization. Yet under this assumption, mod- 
els with an ionization parameter, logC/ ~ —3.5, that 
matches the observed ratios of ions from the same el- 
ement fail to match the column densities of all these ion- 
ization stages unless the metallicity is 5 — 50%, as shown 
in Figure [8[ This indirect measure of the gas metallic- 
ity is affected less by the saturation of H I than might 
be expected. As shown in the right column of Figure [8] 
the column densities of the first and second ions increase 
slowly with respect to iV(H I) through the range where it 
increases from 18.0 to 18.8, because over this range the 
H is becoming increasingly neutral and is added to the 
line integral through the cloud at a much higher rate, per 
pathlength, than species such as Mg II or Fe II for which 
the ionization fractions change little over this range. For 
a given U, these cumulative column density curves shift 
up and down in those panels of the figure in direct pro- 
portion to the assumed cloud metallicity, which has been 
adjusted to best match the greatest number of column 
densities. The best fitting model has O.IZq, but values 
from 0.05 — Q.5Zq still provide decent matches for most 
ions. Even though the metallicities estimated directly 
from the detection of O I are highly uncertain because 
of the saturation of the H I, they are consistent with 
these models at the higher end of the permitted range of 
7V(H I). 

Finally, we have assumed a nominal density of uh = 
0.01 cm"-^, which yields cloud sizes of 1 — 3 kpc. We 
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therefore have a model in which the two strongest com- 
ponents in the LLS absorber are consistent with dif- 
fuse, possibly low-metallicity clouds of total gas column 
logiV ~ 19 — 20 and with kpc scale that are ionized by 
the stars in their nearby galaxies. These absorbers sug- 
gest that star- forming galaxies with /esc = 0.1 may be 
able to photoionize some of their own halo gas under the 
right conditions. 

In light of the appearance of many ionization stages in 
this absorber, e.g. both O I and O VI, we must consider 
that some or all of the gas arises in warm or hot coUi- 
sionally ionized material, such as might arise in a shock 
or conductive interface between hot and cold gas. For 
purely collisionally ionized gas without a photoionizing 
field, we can use individual ionization ratios to constrain 
the temperature. The ratios of the second ions to the first 
ions of C, N, Si, and Fe (e.g. N III to N II) aU prefer a 
temperature of between 20000 — 80000 K. Over this range 
the neutral fraction of H drops from ^ 0.1 to 10"'*, so 
if collisional ionization dominates in these absorbers, the 
total column density of H will exceed logiV(H I) — 19.5 
and could be substantially higher. However, the observed 
ratio of Mg II / Mg I = 1.8 dex prefers a temperature 
T < 20000 K, and is one of the better determined ratios 
in components B and C. Thus, no single temperature in 
CIE describes all the detections. By contrast, the pho- 
toionization model above appears to be able to account 
for these ratios altogether, with the significant exception 
of the O VI. 

An important caveat to this modeling is that we have 
considered ionization mechanisms in isolation, rather 
than possible combinations of them. Such a model might 
explain the co-existence of Mg II and O VI, for instance, 
if the latter is collisionally produced while the former 
arises in a photoionized cloud. Such a model will require 
more parameters than either model on its own, but may 
be tenable. We will discuss the possible "multiphase" 
nature of the O VI absorber below. 

Indeed it is impossible to reconcile O VI and low ions 
such as Mg I into the same GIE or photoionization model, 
and both of these appear in components B and C. Since 
the O VI is not readily reproduced by either of the pho- 
toionization models considered above, it probably does 
arise in collisionally ionized gas. If so, the correspond- 
ing H I may be too weak to separate from the stronger 
H I from the photoionized component. If the O VI ex- 
ists at its peak ionization fraction fovi = 0-2 in CIE 
at T = lO^ ''^ K, the corresponding H I has N(H I) 
~ (2.5/Z) X 10^^ cm~^ where Z is the metallicity rel- 
ative to solar (thus 2.5 x lO^'^ cm^^ at O.IZ0). Only 
if the assumed temperature drops to 100, 000 K, where 
fovi — 10~^, does the associated H I exceed 10^^ cm~^ 
and so become comparable to the observed values. How- 
ever, at this temperature the H neutral fraction is also 
~ 10~^ and the total H column is an implausibly large 
10^^ cm~^. We therefore conclude that a collisional ori- 
gin for the O VI (and possibly other high ions for which 
we do not have coverage) is possible, and that the cor- 
responding H I would go undetected under most circum- 
stances. Such a case might occur if the O VI arises in the 
interface layer between a photoionized cloud and a hot, 
diffuse halo, as has been proposed for the high-velocity 
cloud population of the Milky Way ( Sembach et al.|2003 



Component A 




Fig. 9. — Parameter space for photoionization models corre- 
sponding to Component A. The most informative constraints are 
imposed by the column densities of C II (blue), Mg II (orange), 
S III (black), and S IV (red). 



Fox et al. 2004 1. As these interfaces lie at intermediate 
temperatures between hot and cold material, they may 
account for some portion of the intermediate ionization 
stages as well. Were we able to account for this properly, 
it might change our conclusion about the photoioniza- 
tion models described above. However, these interface 
models typically produce less O VI than is observed, as 
discussed more below. 

We are therefore forced to conclude that these compo- 
nents consist of complex multiphase material, and that 
they may arise from some combination of photo- and 
coUisional-ionization. The essential point of these ion- 
ization modeling exercises is not that a particular model 
gives a perfect fit to the observed ion ratios; the uncertain 
N(H I), component blending, and line saturation ensure 
that no perfectly robust modeling is tenable. However, 
ionization models with a range of plausible conditions 
all point in the same direction; that the bulk of the gas 
traced by the H I and the metal-line absorption is highly 
ionized, with a total H column density exceeding that 
seen in H I by a factor of 10 - 100. There is still more 
highly ionized gas, traced by O VI in all the components, 
that may indicate interfaces between this photoionized 
halo gas and its environment. While we cannot com- 
pletely rule out collisional ionization for the bulk of the 
strong components B and C, a photoionization model 
incorporating UV light from star forming regions in the 
nearby galaxies appears plausible if the cloud has density 
uh — 10~^ cm~^, lies 20 - 50 kpc from the associated 
galaxies, and if those galaxies propagate ~ 10% of their 
ionizing photons into the halo. 

5.2. Component A 

This component exhibits a wide range of ionization 
stages, from Mg II to O VI. This range of species can- 
not be accommodated within a single photoionization or 
collisional ionization model, for the same reasons that 
applied to B and G above; it is clearly "multiphase". 
Following the same procedure as above, we attempted to 
fit as many species as possible into a single model to esti- 
mate what properties the cloud might have under various 
physical scenarios. 

We first ask whether the combination of a galactic star- 
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forming spectrum with the diffuse extragalactic back- 
ground will provide a suitable model, as it did for com- 
ponents B and C. Since we have a robust measurement 
for logiV(H I) = 16.5 in this case, we can restrict our 
attention to models that give that value. These results 
appear in Figure [9j Most regions of this parameter space 
with log U > —3 (from the combination of the HM and 
S99 spectra) satisfy the lower limits on C III and Si III. 
The measured column densities of C II, Mg II, S III, and 
S IV are constraining, and mutually consistent in a small 
region centered around the filled circle at \ogU — —2, 
Z = 0.2 Zq. This model has a total H column of 
log A^(H) ~ 19, and thus a neutral fraction of only 10~^-^. 
However, the ionization parameter is 10 - 30 times higher 
than for B and C, suggesting that this component has 
some combination of lower density and/or proximity to 
the ionizing source. We therefore conclude that this pho- 
toionizing model is plausible for all the species except the 
O VI. 

The O VI detection suggests that some of this gas may 
reside in a hotter, coUisionally ionized component. We 
pursued this hypothesis by adjusting the temperature 
and metallicity of a purely coUisionally ionized model 
cloud to achieve matches to the observed ratios and col- 
umn densities. No single temperature is clearly indi- 
cated; the well-measured ratio of S IV to S III is matched 
by a model with 80000 K, but this model greatly under- 
produces both Mg II and O VI, which has its highest 
ionization fraction at ~ 300000 K. Single-temperature 
models that give O VI column densities above 10^'* cm~^ 
typically give Mg II column densities 100 x lower than we 
measure here. A model with logT = 5.3 gives a ratio of 
O VI to S IV that is near the observed value, but overpre- 
dicts the observed ratio of S IV to S III. Where therefore 
conclude that no single-temperature collisional ionization 
model can explain the absorption by component A; dis- 
tinct models are needed for the low ionization gas up to 
C III, Si III, and S III, while a hotter, coUisionally ionized 
phase could account for the O VI and possibly the S IV. 
In a combined model, the neutral H fraction associated 
with the collisional component is likely quite low with 
respect to the measured log A^(H I) = 16.5, which would 
trace predominantly the photoionized component. We 
thus have no constraint on the metallicity of the hotter 
component. 

Thus we find that Component A can also be accommo- 
dated into a simple scenario in which low-metallicity gas 
in the halo(s) of the nearby galaxies is ionized by their 
star formation. Such a simple scenario can account for 
the observed column densities with the significant excep- 
tion of the O VI, for which a coUisionally ionized phase 
may be preferred. 

5.3. Component D 

Component D has a column density logiV(H I) = 14.8, 
and exhibits strong absorption in O VI, saturated C III 
and Si III, and the problematic transition of N III. Signif- 
icantly, no Fe III is seen here as in B and C. With fewer 
detected ions, models are more loosely constrained. In 
this case, a single temperature CIE model can account 
for the column densities of C III, Si III, N III, and O VI 
in a narrow temperature range around logT ~ 5.25, pro- 
vided the total H column density is logAf(H) ~ 20 and 
the metallicity is solar. However, this model conflicts 
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Fig. 10. — Parameter space for photoionization models corre- 
sponding to Component D. The strong detections of C III and 
Si III exclude models to the left of the green and blue lines, while 
the upper limit on Mg II requires models below the orange line and 
the upper limit on S IV requires models to the left of the red line. 

with the observed doppler parameter 6 = 16 km s~^ ob- 
tained by consistent fits to the profiles of Ly/3 and Lyj. 
At logT = 5.3, the doppler parameter for H should be 
no less than b = ^j2kT /mu — 57 km s~^, but the line 
profiles prefer a value T < 20000 K, which corresponds 
to 6 = 16 km s~^. We therefore rule out a single CIE 
model for all the detected absorption in this component. 

If we consider a photoionization model for this compo- 
nent, we find that the number of detected ions is gener- 
ally insufficient to constrain the ionization parameter to 
the degree that was possible for A, B, and C. Any model 
with \ogU > —2.5 gives the correct column densities of 
C III and Si III, but tends to produce too much Mg II, 
which is undetected, unless the metallicity is less than 
0.3 sol ar. This region of parameter space is shown in 
Figure [To] The allowed area with logJ7 = 2.5 to —1 and 
below 25 — 40% solar produces enough C IV and Si III 
without producing too much Mg II or S IV. The total gas 
column densities N(H) range from 10^® — 10^^'^ cm~^ as 
\ogU increases from —2.5 to —1. This component would 
appear to have a very large ionization correction, with 
only 10^'^ or lO"** of the hydrogen in H I. As for the other 
components analyzed above, we find that the combined 
ionizing radiation can explain the detected absorption, 
but only if the gas is highly ionized and has a subsolar 
metallicity. 

5.4. What Creates the O VI? 

None of the analysis just presented answers convinc- 
ingly the question of what creates the highly ionized gas 
traced by the O VI. Models with purely photoionized gas 
do not produce enough O VI, or if they do they do not 
also explain the lower-ionization gas. It is the coexis- 
tence of this strong O VI, together with the low ions, 
that make this absorber "multiphase" . 

We must consider the possibility that the highly ion- 
ized gas in this system - the O VI and perhaps some or 
all of the S IV, S III, and N III - arise in non-equilibrium 
situations such as shocks, conductive interfaces between 
hot and cold material, or simply cooling ga s which is 
out of ther mal and/or ionization equilibrium. Fox et al. 
( |2004| and |Indebetouw fc ShuU] (2004a,b^ have compiled 
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models for many of these phenomena, and mapped out 
their predicted column-density ratios in Si IV, C IV, N V, 
and O VI. Unfortunately our COS spectra have coverage 
of only the latter two ions, and N V is not detected. 
The ratio N(N V) / N(0 VI) provides only an upper 
limit for the four component groups, which ranges from 
< —0.6 to < —1.0. These limits provide no constraint, 
since they are consistent with shock ionization, radiative 
cooling, conductive interfaces, and tu rbulent mixing lay- 



ers, as compiled by Fox et al. (2004). The uncertainty 



of comparing models to data in this fashion is compli 
cated by the possibility that the gas under study has 
non-solar relative abundances of the diagnostic elements 
C, Si, N, and O, which allows the regions occupied by 
these various mode ls to move within the column density 
ratio space. While Fox et ah] (2004) are able to correct 
these abundance ratios based on independent measure- 
ments of relative abundances, we do not have that much 
constraint here. So we are unable to produce quanti- 
tative constraints on ionization mechanisms using this 
technique. 

This absorber is notable not only for its very high total 
column density on O VI, but also for how it produces this 
high value. We do not see a monolithic, saturated com- 
ponent in O VI, thermally broadened into a single profile 
in hot gas, nor do we see O VI associated only with the 
strongest or weakest H I. Rather, we see a complex profile 
in which the O VI absorption is spread over nearly 400 
km s~^ and breaks into at least two and possibly more 
distinct components, which line up with components at 
similar velocities in lower ions. With a sightline that ex- 
tends over > 100 kpc through the halos of two galaxies, 
it seems natural to guess that the O VI arises in multiple 
distinct objects within the halo which may or may not 
be associated with the individual galaxies. Strong shocks 
(v ~ 600 — 2000 km s~^) c an exhibit column d ensities in 
the range of our detections ( Gnat fc Sternbergj2009,) , but 
imply post-shock temperatures of T ~ 5 x 10" K and H I 
thermal linewidths b > 70 km s~^. The strength and 
complex component structure of the H I lines precludes 
a search for such a broad component, which would only 
be detectable in the Lya profile in this Lyman-limit sys- 
tem. Thus while we cannot rule out shocks for some of 
this O VI, we have no strong indication for them either. 
We note that some O VI absorbers reported in the lit- 
erature have corresponding broad H I components that 
could indicate the expected temperature in shocked gas 



( Tripp et al. 2001 Richt er et al. 
^OlOl Savage et a/.^Oll) 



2004 Narayanan et al. 



One possible explanation for the observed kinematic 
and ionization pattern is that the O VI arises in tran- 
sition layers between hot and cold gas, such as might 
occur for neutral or photoionized clumps falling through 
a Galactic halo and interacting with the hot halo gas left 
there by the formation of the galaxy. O VI is believed 
to trace interfaces between stable cold (T ~ 10** K) and 
hot (T > 10^ K) phases of the ISM and circumgalac- 
tic medium because it achieves its maximum ionization 
fraction (T =100-300,000 K) near the peak of the radia- 
tive cooling curve and so is short-lived. In this simple 
scenario, interfaces arise where cold clouds interact with 
the hot medium and are generally unrelated to, or weakly 
dependent upon, the size and/or quantity of gas in the 



cold clumps. Most simply there is just one interface per 
cloud, however large the cloud. 

Quantitative models of conductive interfaces and tur- 
bulent mixing layers typically produce column den fsities 
of O VI in the range ipi^-is ^jj^-^^ lower than we detect 



here (cf . i Bo rkowski et_al 1990[ Slavin et al. 1993 Indebe- 
touw fc Slui ll 2004a|b[ IGnatet al |20ld| , and so multiple 
interfaces must be invoked. 'I'he models of lGnat et al.l 



(20101 find that O VI column densities typically have 
N(0 VI) < 10^^ cm^^ for cold clouds bounded by con- 
ductive interfaces within a hot corona, so that ^ 100 
interfaces are implied by the total column density of the 
J1009-I-0713 absorption-line system. A very large num- 
ber of low-column density O VI components probably 
could be accommodated by our data, since blending and 
thermal broadening would likely make them difficult to 
distinguish from a smaller number of higher column den- 
sity components, which is which is similar to what we 
observe. 

However, the interface scenario implies a simple ob- 
servational pattern: that O VI should appear whenever 
low ions trace colder clumps in the halo, and that the 
strength of O VI absorption near galaxies should increase 
with the number of detected components in the cooler 
material, and correlate with them in velocity space. The 
overall pattern observed so far is that is that the col- 
umn density of O VI is relatively insensitive to the total 
A^(H I) in surveys of IGM absorbers, though the pat- 
terns depend sensitively on the number of defined compo- 
nents, which vary from study to study especially for the 
strongest absorbers. The column density of O VI varies 
in our case by only about 0.2 dex across the four detected 
component groups, while the H I column density varies 
by at least three and possibly four orders of magnitude. 
This striking lack of correlation suggests that the pres- 
ence and total quantity of O VI is not related directly to 
the total quantity of gas. Indeed, this is a generic pattern 
of O VI absorbers that has been noticed in systematic 
surveys fPanfo rth fc Shun|2008[ [Tripp et al.|2008 



Thom 



Ifc Chen 2008). In the cornpilation of U VI absorbers by 
Thom & ChenJ ( |2008[ ), N{0 VI) varies by 1.5 dex over 
5 decades of variation in 7V(H I). This is partly due to 
the dramatically multiphase character of some of the ab- 
sorbers in the Thom & Chen sample, but even the "sim- 
ple" O VI absorbers (systems with well-aligned O VI and 
H I components and simple component structur e with 



no clear evi dence of multiple phases) identified by Tripp 
et al. (|2008 ) show this effect, with N(H I) varying by 3 or 



ders of magnitude while N(0 VI) only changes by 1 dex. 
Taken as a system, the J1009-I-0713 LLS occupies the 
extreme upper right of this diagram, with log iV(0 VI) 
= 15.0 at logiV(H I) ^ 18.5, but taken by components 
the four groupings sc atter over the full range found by 



Thom fc Chenl ( |2008[ ). Thus the comparison of O VI to 
H I is roughly consistent with the interface scenario. 

The low ion Mg II is another cold gas tracer that could 
provide information on the interface scenario. We have 
successfully built a photoionization scenario that works 
well for the detected Mg II and other low ions. The sim- 
plest interface scenario suggests that each "cold cloud" 
could contribute two interfaces to the sightline. We have 
grouped the major kinematic components in this ab- 
sorber into four ranges, but if components are generously 
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TABLE 3 
Literature on OVI-bearing LLSs 



# 



Siehtline 



^abs 



1 

2 
3 
4 
5 
6 
7 
8 
9 
10 
11 



HE0153-4520 
PKS0312-77 
PKS0405-123 
PG1116+215 
PG1216+069 



3C 351.0 
PHL1811 

J1009+0713 



0.022601 
0.2026 
0.16692 
0.13847 
0.00632 
0.12360 
0.28189 
0.22111 
0.07765 
0.08093 
0.3558 



i6.6ii«:- 

18.22l«-i«3 
16.45 ±0.05 
16.20 ± 0.04 
19.32 ±0.03 

> 15.95 
16.70 ± 0.04 

> 17.0 
16.03 ±0.07 
17.98 ±0.05 
18.0 - 18.8 



14.21 ± 0.01 
14.95 ±0.05 
14.72 ±0.02 
13.68l°:J° 

< 14.26 
14.83 ±0.10 
14.02 ±0.02 
14.27 ± 0.04 
13.56 ±0.10 

< 13.59 
15.0 ±0.2 



1 
2 

3-7 

8 

9 
10 
10 
10 
10 
11 
12 



)0^; (3) pieiTfe 
: (6)|Lehner et aL 



References. — (1) Savage ct al. (20111; (2) Lclmcr ct al. (20091, 
[Prochaska (2000l); (4) fProchaska ct al. (20041; (5) Willigcr ct al.^ (,2006|rre) [ 
^PJ7); (7) Savage ct al. (201U); (8) Scmbach ct al. (20U4); (gj jl'ripp et al. | ( |2UU5"t ; ( lU) 
iTTipp ct al. ( 2008 ); (ll) _Jcnkins ct al.^ ( ,2U05) . (12) this paper. 



defined to maximize their number for the Mg II there 
could be approximately nine components over groups A 
- C, but none for group D. Thus at most 20 interfaces are 
expected in the O VI, not the > 100 that are implied by 
the observed column density and the interface models. 
While there could be numerous cold clouds and twice as 
many interfaces that fall below even the stringent detec- 
tion limits of our HIRES data, there would need to be 
a large number of undetected Mg II to provide enough 
interfaces to reproduce the O VI. 

The features of this system recall the properties of 
some other low-redshift, O Vl-bearing LLS that have 
been intensively studied to which we can compare and 
contrast our results. These systems are tabulated in 
Table 3, where we list all the systems in the literature 
with well-studied UV spectra covering O VI from HST 
or FUSE and having logA^(H I) > 16 to within obser- 
vational errors. The listed number of components is the 
number of H I and/or metal-line components given in the 
references cited. Here we have assigned 9 components to 
the J1009+0713 LLS based on the number of observed 
Mg II components in the HIRES data, which has compa- 
rable resolution to the STIS/E140M data used in most 
of the other cases from the literature. This compilation 
of results suggests two possible patterns of relevance to 
the J1009-I-0713 LLS and the origins of strong O VI ab- 
sorbers. First, the total quantity of O VI is not well 
correlated with the column density of H I. The former 
spans less than 1 order of magnitude while the latter 
covers more than three dex. This poor correlation is one 
of the central p uzzles in the origins of the O VI absorbers 



(cf. Figure 6 of Thom & Chen 
system column density ot U v' 



( |2008[ )). Second, the total 
appears to correlate with 



the number of detected components in the system, such 
that > 4 components gives \ogN{0 VI) > 14.5, while 
systems with a smaller number of components have lo wer 
column density. This effect is clearly seen in Figure |1H 
where we plot the O VI column density over the number 
of components from Table 3. 

While the sample is still small, there is a marked corre- 
lation that is well outside the typical error on individual 
measurements (conservatively assigned as ±0.1 dex and 
shown at lower right). We apply a Spearman's rank cor- 
relation test and find that we can reject at 99.7% con- 
fidence the null hypothesis that there is no correlation 



between the O VI column density for these systems and 
the number of kinematic components they exhibilr^ We 
emphasize that this is not a homogeneously selectea sam- 
ple, though we have tried to gather all reported systems 
that meet the selection criteria. Nevertheless we find at 
least a strong suggestion that a large number of com- 
ponents leads to a strong O VI absorber. This correla- 
tion could indicate that the O VI is generally associated 
with the low-ionization material even if it does not re- 
side in the same physical layers of the gas. This general 
behavior is expected in the interface scenario discussed 
above, though the problem of how just a few interfaces 
give log A^(0 VI) > 14.5 remains. 

The J1009+0713 system and other systems like it thus 
offer mild but not conclusive evidence in favor of the in- 
terface scenario. The broadening of strong O VI into 
many detectable components and a correlation with a 
large number of velocity components in the H I and low 
ions are consistent with the basic expectations of the in- 
terface scenario. But the total quantity of O VI substan- 
tially exceeds the expectations of the interface models, 
given the number of apparent cold gas clouds along the 
sightline. Three possible solutions to this puzzle are that 
a large number of cold clouds go undetected but their 
interfaces are seen, or the O VI from a single interface 
is well above that calculated in models, or the interface 
model does not apply. It is still possible, perhaps even 
likely, that some fraction of the O VI we observe is pro- 
duced in interfaces, along with additional contributions 
from another mechanism. This system shows that fur- 
ther measurements of O VI and Mg II can provide more 
information about the correlation of hot gas with low- 
ionization gas for the same sample of absorbers. Such a 
study is planned as part of our larger survey. 

6. SUMMARY AND DISCUSSION 

We have examined the ionization, metallicity, and 
association with galaxies of a newly-discovered strong 
intervening O VI absorption-line system discovered in 
our COS data on J1009+0713. This system exhibits 
logA^(0 VI) = 15.0 spread kinematically over 400 km 

The correlation is still strong if we assign only the 4 compo- 
nents detected at COS resolution to the J1009+0713 LLS; the data 
then reject the null hypothesis at 99.3% confidence. 
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Fig. 11. — A correlation of the total system column density of 
O VI versus the number of reported kinematic components for the 
literature compilation reported in Table 3. The typical error bar is 
generously assigned to be ±0.1 dex and is displayed at lower right; 
the apparent trend is well outside this margin of error. 

s^^. It appears to be associated with at least two galax- 
ies at projected separations of 14 and 46 kpc from the 
sightline that have redshifts coincident with the detected 
gas. The system includes two similar LLS-strength com- 
ponents separated by 60-80 km s^^ with two weaker asso- 
ciated components 60 - 100 km s~^ away from these. The 
direct line measurements only constrain the gas metal- 
licity in the two strongest components to Z = 0.1 — IZq 
owing to saturation of the H I. However, the detected 
ion ratios are well-matched by a photoionization model 
in which simple clouds of total column density log Af(H) 
~ 20 and 5 — 50% solar metallicity are ionized to a 1 - 
10% neutral fraction by a radiation field from star form- 
ing regions at a distance of 20 - 100 kpc. In short, it 
appears that the bulk of this absorber traces gas that re- 
sides in the common halo environment of these galaxies 
and is ionized by their ongoing star formation. Simi- 
lar modeling applied to the two outlying components, 
with logiV(H I) — 14.8 and 16.5, are also consistent with 
their detected absorption for models with \ogU ~ — 2 
and ~ 25% solar metallicity. These outlying clouds also 
appear to be relatively highly ionized, with neutral frac- 
tions of 0.1 - 1% even in the gas traced by low ions, but all 
with an additional, difficult-to-constrain, highly ionized 
component traced by O VI. 

This system exhibits two galaxies ranging over two 
magnitudes in the rest-frame optical, the brightest 
(170_9) has one with M* = 6 x lO^M© in start[^ Galaxy 
86_4 is 2 mag fainter with a total star formation rate of 
only about 1/10 that of 170_9. Based on emission-line 
measurements and its inferred stellar mass, the brighter 
galaxy 170_9 appears to be consistent with solar metallic- 
ity, while 86_4 is either consistent with solar metallicity 
or with 25% solar (depending on which branch of the 
R23 indicator is adopted) . Since none of the detected gas 
components appear to require a solar metallicity (though 
B and C are consistent with solar if the photoionization 

The galaxy 92_7 appears in the SDSS photometric database 
with r = 22, and appears distinctly in the WFC3 image, but we 
have not yet obtained its redshift. It may also be associated with 
this absorption-line system. 



modeling is discounted), we do not have positive evi- 
dence that any of the detected gas was brought into the 
sightline by outflows from the galaxies. 

Several features of specific systems are worth compar- 
ing to the J1009-f 071 3 LLS. The system at z = 0.167 to- 
ward PKS0405-1219 ( |Chen fc Prochaska,|2000| hes ~ 70 
kpc from two massive galaxies, is at least 0.1 solar and 
likely higher metallicity, and exhibits strong absorption 
from high ions up to N V and O VI at log iV( VI) = 
14.6. The z = 0.0809 LLS toward PHL 1811 (jJenkins 



et al.|2005 ) lies 34 and 87 kpc from two nearby L"^ galax- 



les, and has nearly a solar metallicity, but is not asso- 
ciated with highly ionized gas. Both of these absorbers 
show evidence of a predominantly ionized gas, though 
the PHL 1811 system also exhibits absorption by low 
ions. While our simple photoionization models shown 
above indicate a sub-solar metallicity for the J1009+0713 
system, this model-dependent value may not reflect the 
true gas metallicity in the LLS components: the O I 
measurements are consistent with solar metallicity. The 
PKS0405-1219 system is a particularly interesting com- 
parison to the J1009+0713 system, since the O VI in 
both cases is stronger and more kinematically complex 
than would be expected from the lower ionization species. 
While is it difficult to draw general conclusions from only 
three cases, it appears that galaxy halos may commonly 
host LLS absorption, that it can be quite enriched in 
metals to nearly solar abundances, and that it can be 
associated with large quantities of highly ionized gas. 
Based on the PKS0405-1219 and J1009+0713 systems, 
we speculate that LLS-strength absorption may be neces- 
sary to produce the strongest intervening O VI absorbers, 
logiV(0 VI) > 14.5. The general incidence of LLS in 
galaxy halos and its relationship to hot gas will be ad- 
dressed by our larger survey, which should also be able to 
assess the gas metallicity in many cases. High metallicity 
in dense, cool gas located at > 50 kpc from the nearest 
galaxies may indicate a major role for tidal stripping of 
dwarf satellites, or robust outflows from large galaxies or 
their dwarf satellites. 

This sightline passes by a galaxy with nearly and 
a dwarf galaxy possibly in its halo or interacting with 
it, which may have low metallicity. Their total velocity 
separation is only 25 km s~^, which suggests they are 
gravitationally associated. It is tempting to draw com- 
parisons of this system to the Milky Way and its Mag- 
ellanic Clouds. If the smaller galaxy has the metallic- 
ity of 25% indicated by their lower-branch R23 measure- 
ments, it could be associated with the A and/or D com- 
ponents in an outflow or tidal-stripping scenario. If so, 
some of the detected absorption components may arise in 
gas stripped from this dwarf galaxies, like the Magellanic 
Stream. 

Alternatively, all the detected absorption components 
could have sub-solar metallicity and could trace low- 
metallicity, infalling material not unlike the large HVC 
complexes of the Milky Way. This inference and the high 
ionization state of the detected gas would imply that in- 
falling, metal-poor gas enters halos and/or galaxies in 
ionized form, and is possibly associated with still more 
highly ionized material traced by O VI that arises ei- 
ther in the hot diffuse galaxy halo itself or in interface 
zones of intermediate temperature between that hot am- 
bient medium and the cooler infalling clouds. This see- 
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nario has been well developed for the Milky Way - further 
tests of its validity for external galaxies requires better 
knowledge of gas metallicities, kinematics in association 
to galaxies, and thus more well-studied cases. 

In our combined COS and HIRES data we have used 
two important diagnostic lines for gas in galaxy ha- 
los - both O VI and Mg II - that are typically not 
available in combination for the practical reason that 
they require both U V and optical data and galax ies at 
z > 0.1. Recently, iBarton & Cookel (|2009|) and ^heK 



et al. (2010) have systematically addressed the number 
density ot Mg II absorbers in the ~ 100 kpc regions sur- 
rounding galaxies using galaxies and QSO sightlines se- 
lected from SDSS. They find that the equivalent width 
of Mg II increases at low impact parameter in a well- 
defined trend that has a detectable dependence on galaxy 
luminosity for those galaxie s that exhibit Mg II ab sorp- 
tion. However, both ^Bowen fc Chelouche (2011) and 
Gauthier et al. (2010) have found lower covering frac- 
tions of strong Mg II surrounding more massive galaxies 
(LRGs) that indicate the cool gas covering fraction de- 
pends, perhaps strongly, on galaxy type and/or environ- 
ment. For comparison with these samples, we obtained 
a total = 1250 ± 15 mA for the full LLS analyzed 

here, including 188 mA in component group A, 604 mA 
in component group B, and 472 mA in component group 
C. These findings are all within the range of the distri- 
butions found by these other studies. However, since we 
also detect H I, O VI and other multiphase UV ions in 
these component groups, we can draw additional impli- 
cations that are not possible based on Mg II alone. 

First, we have used the detections of Mg II and the 
ionization stages of C, N, O, Mg, S, Si, and Fe to infer 
that the absorbing clouds have only a small portion, 0.1 
- 1% of their gas in the neutral phase. This finding im- 
plies that the gas mas s traced by t hese a bsorbers could 
be significant. Indeed, Chen et al. (2010) inferred Mg II 
cloud masses of ^ 2 x 10* Mq and a total baryonic mass 
in halos of 3 x 10^ M©, using the important assumptions 
of 10% solar metallicity and a 10% ionization fraction 
for Mg II. These values are within the range permitted 
for the J1009-1-0713 clouds, where the ionization correc- 
tions and metallicities are better constrained by multiple 
ionization stages. Thus our findings indicate that the as- 
sumption of a significant mass correction to the observed 
Mg II is reasonable, and such absorbers could harbor sig- 
nificant mass. Moreover, the O VI suggests a possibly 
significant contribution of mass over and above that pro- 
vided by Mg II systems. Finally, it may prove important 
to measure halo absorbers in both O VI and Mg II; com- 
ponents A, B, and C here hint that O VI and Mg II can 
have some relationship but may vary in their ratio by over 
an order of magnitude. Our larger survey was designed 
with the goal of detecting both lines in a significant sam- 
ple of galaxies over a ran ge of mass. We will thus b e able 
to tes t the findings of [Barton fc Cooke| (|2009) and Chen 



et al. (2010j) that Mg II cor relates with galaxy lu mmos 



ity, and also the find i ng by Bowen & Chelouche (20111 
and Gauthier et al. (2010) that massive galaxies have 
less Mg II, while also exploring these relations in O VI 
for the same galaxies. These previous results, considered 
together with hints of a Mg II / O VI relationship for the 
J1009-I-0713 system, suggests this will be a fruitful line 



of research. 

While it is among the strongest known O VI absorbers, 
this system at least superficially resembles the highly- 
ionized, multiphase, O Vl-be aring HVCs surro unding 



the Millky Way. For example, |Fox et aE] (p005| report 

ra 



conditions very similar to those we find here for two 
HVC complexes toward HE0226-4110 and PG0953-H414 
in the Milky Way halo. Specifically, their ionization 
modeling finds logJ7 ~ —3.5, subsolar metallicity, and 
lognn =i -2, but for clouds with logiV(H I) - 16-16.5. 
Like them, we have found consistent photoionization 
models for the present system that yet do not explain the 
observed abundance of O VI. Our findings also resonate 
with the recent survey of ionized si l icon i n th e Milky 
Way HVC p opulation by |Shull et al] ( |2009[ ) and [Collins 
et al. ( 2009 1 , who find typical values of log U ~ —3 and 



10 - 30% solar metallicity or lower in these clouds with 
neutral fractions of ~ 1%. These favorable comparisons 
suggests that highly-ionized, possibly metal-poor gas re- 
sides in the halos of star forming galaxies as it does in 
the Milky Way - a key motivation for the survey that 
discovered this absorber. A more speculative but inter- 
esting inference can be drawn from the column density 
and ionization of the main components B and C. Were 
they exposed to a lower radiation field, because of a lower 
galactic f^sc or SFR, they could instead appear as a sub- 
DLA or DLA with iV(H I) 10^*' cm'^ located at least 
20 kpc from the nearest galaxy, and might have total 
column densities similar to the classical MW HVC com- 
plexes. We note that recent surv eys for 21-cm absorption 



affiliated with near by galaxies (Borthakur et al. 2010 
Gupta et al. 20101 have found 21-cm absorbers to be 



quite rare m galaxy halos. While this could be an in- 
dication that sub-DLAs and DL As are rarely found at 



such large inipact parameters (cf. Meiring et al. 2011b 



Peroux et al.[[2011 1, this remains an open question since 



21-cm absorption depends on spin temperature, and sub- 
DLAs and DLAs might be missed in 21-cm absorption 
surveys if the gas is too warm. Our larger survey will 
address this problem by assessing the incidence of LLS, 
sub-DLAs, and DLAs in galaxy halos. 

While the similarity of this absorber to some well- 
studied Milky Way HVCs provides an encouraging link 
between the HVCs and galaxy halo gas at large, this 
system also illustrates some important limitations that 
must be overcome to fully map out diffuse gas in ha- 
los. First, H I column densities in the range that gives 
an LLS but not a DLA, as here, are extremely difficult 
to constrain to a narrow range. This is especially true if 
the gas obviously occupies several clearly separated kine- 
matic components that are blended in the Lyman series, 
as shown by our HIRES data. This unfortunate circum- 
stance makes metallicities hard to estimate, even in cases 
where conditions are otherwise favorable for obtaining 
metal-line column densities, as they are here. While LLS 
are rare in the IGM at large, they may be quite common 
close to galaxies where total gas densities are high, even if 
the gas is typically highly ionized. This system also illus- 
trates vividly that relating halo gas to its host galaxies is 
complicated by projection effects and by multiple galax- 
ies lying along the sightline. We clearly cannot associate 
the absorption, especially in the O VI, with any of the 
three galaxy candidates in particular. 

Finally, this system illustrates that detection of many 
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ionization stages can constrain the physical parameters 
of the gas, but only for well-defined hypothetical phys- 
ical scenarios, e.g. pure photoionization or pure col- 
lisional ionization. Combinations of such models, or 
non-standard scenarios such as interfaces and/or non- 
equilibrium calculations rapidly become too complex to 
test easily. Even for simple scenarios it is sometimes 
difficult to identify unique explanations for detected ab- 
sorption. These considerations suggest that a large, sys- 
tematic survey to assess possible variations could help by 
smoothing out variations along individual lines of sight. 
Our systematic survey of gas in galaxy halos, which 
serendipitously discovered this absorber, aims to estab- 
lish a set of empirical facts about the content of gaseous 
halos. Yet multiphase, kinematically complex absorbers 
like this one pose a stiff challenge to our ability to build 
theoretical models that are up to the task of relating 
complex physical processes in galaxy halos to observable 
quantities. Progress in this area could significantly ad- 
vance the cause of understanding galaxy formation at 
large. 
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